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Balance functions have been measured for charged particle pairs, identified charged pion pairs, 
and identified charged kaon pairs in Au+Au, d+Au, and p+p collisions at y^SNN = 200 GeV at the 
Relativistic Heavy Ion Collider using the STAR detector. These balance functions are presented 
in terms of relative pseudorapidity, A?7, relative rapidity. Ay, relative azimuthal angle, A(/), and 
invariant relative momentum, ^inv In addition, balance functions are shown in terms of the three 
components of qinv' ^long, ^out, and ^side- For charged particle pairs, the width of the balance 
function in terms of Ar] scales smoothly with the number of participating nucleons, while HIJING 



and UrQMD model calculations show no dependence on centrality or system size. For charged 
particle and charged pion pairs, the balance functions widths in terms of Ar] and A^ are narrower in 
central Au+Au collisions than in peripheral collisions. The width for central collisions is consistent 
with thermal blast-wave models where the balancing charges are highly correlated in coordinate 
space at breakup. This strong correlation might be explained either by delayed hadronization or 
by limited diffusion during the reaction. Furthermore, the narrowing trend is consistent with the 
lower kinetic temperatures inherent to more central collisions. In contrast, the width of the balance 
function for charged kaon pairs in terms of Ay shows little centrality dependence, which may signal 
a different production mechanism for kaons. The widths of the balance functions for charged pions 
and kaons in terms of ^inv narrow in central collisions compared to peripheral collisions, which may 
be driven by the change in the kinetic temperature. 

PACS numbers: 25.75.Gz 



I. INTRODUCTION 

The study of correlations and fluctuations can pro- 
vide evidence for the production of a strongly interacting 
quark-gluon plasma (QGP) in relativistic heavy- ion colli- 
sions [imQl] . Various theories predict that the production 
of a QGP phase in relativistic heavy-ion collisions could 
produce significant event-by-event correlations and fluc- 
tuations in temperature, transverse momentum, multi- 
plicity, and conserved quantities such as net charge. 

One such observable, the balance function, maybe sen- 
sitive to the correlation of balancing charges [27]. For 
instance, for every particle of momentum p, there must 
be an anti-particle of momentum p' with the opposite 
charge. By means of a like-sign subtraction, the bal- 
ance function can produce the distribution of relative 
momentum, q = p — p' ^ between the balancing charges. 
Balance functions are sensitive to the mechanisms of 
charge formation and the subsequent relative diffusion 
of the balancing charges [27] . Balance functions are also 
affected by the freeze-out temperature and radial flow 
[23 • Remarkably, balance functions for central collisions 
have been shown to be consistent with blast- wave models 
where the balancing charges are required to come from 
regions with identical collective flow [30]. The inferred 
high degree of correlation in coordinate space has been 
postulated as a signal for delayed hadronization |27j , 
which would not allow charges the opportunity to sep- 
arate in coordinate space. The idea is that in central 
colHsions a deconfined system of quarks and gluon is cre- 
ated, which cools and expands. The observed balanc- 
ing charges are then created when the deconfined sys- 
tem hadronizes, which reduces the effects of expansion 
and diffusion on the correlation of the balancing charges. 
The same arguments were used in discussing charge fluc- 
tuations [5]. Additionally, the same correlations would 
ensue if the charges were created early (on the order of 
1 fm/c), but due to very limited diffusion, remained cor- 
related at breakup. Thus a narrowing of the balance 
function in central collisions implies delayed hadroniza- 
tion. We have previously presented results for balance 
functions from Au+Au collisions at -s/snn = 130 GeV 
for all charged particles and for identified charged pions 
|23| . We observed that the balance function narrows in 
central Au+Au collisions for all charged particles and for 



identified charged pions. 

UrQMD (Ultra-relativistic Quantum Molecular Dy- 
namics, version 2.3) [31] is an example of a model where 
charges are created early and there is significant diffu- 
sion in the subsequent evolution of the system. Indeed, 
balance functions in terms of relative pseudorapidity or 
relative rapidity predicted by UrQMD do not exhibit nar- 
rowing in central collisions (see Section V). Other mod- 
els have been applied to predict balance functions. One 
model is based on a blast- wave and includes a thermal 
model with resonance decay [32|- This model cannot ex- 
plain the narrowing of the balance function in central 
Au+Au collisions at -^^snn = 130 GeV. Another model 
attributes the narrow balance functions observed for cen- 
tral Au+Au collisions at ^/sNN = 130 GeV to quark- 
antiquark coalescence |33| . 

Recently, the system size and centrality dependence of 
the balance function for all charged particles has been 
studied at ^/sNN = 17.3 GeV for p+p, C+C, Si+Si, 
and Pb+Pb collisions [3^. The balance function for all 
charged particles narrows in central Pb+Pb collisions at 
17.3 GeV and the widths of the balance functions for 
p+p, C+C, Si+Si, and Pb+Pb collisions scale with the 
number of participating nucleons. The rapidity depen- 
dence and incident energy dependence of the balance 
function for all charged particles have been studied for 
Pb+Pb collisions from ^/sNN = 6.3 GeV to ^/sNN = 
17.3 GeV in Ref. [35]. The balance function is ob- 
served to narrow in central collisions for midrapidity, but 
does not narrow at forward rapidity. The authors of Ref. 
[35] show that the narrowing of the balance function in 
terms of A?^ in central collisions can be explained with 
the AMPT (A MultiPhase Transport) model incorporat- 
ing delayed hadronization, while models such as HIJING 
and UrQMD fail to reproduce the observed narrowing. 
We have recently presented a study of the longitudinal 
scaling of the balance function in Au+Au collisions at 
./I^ = 200 GeV [36]. 

In this paper, we present new results for the balance 
function from p+p, d+Au, and Au+Au collisions at 
^/sNN = 200 GeV. These results have significantly bet- 
ter statistical accuracy than our previous measurements 
for Au+Au collisions at s/snn = 130 GeV and define 
the system size dependence of the balance function at 
^/sNN = 200 GeV. We present balance functions for all 



charged particles, charged pions, and charged kaons. We 
also show the balance function in terms of several dif- 
ferent variables that each have different sensitivities to 
different physical effects. We compare our results with 
current theoretical predictions. 

The balance function is calculated as: 
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where A^ in the case of identified charged pion pairs 

denotes the density (number divided by bin width) of 
identified charged pion pairs in a given range, e.g. rel- 
ative rapidity Ay = |^(7r+) — ?/(7r~)|, and similarly for 

A++, A__, and A ^. The terms A^ , A++, A , and 

A_+ are calculated using pairs from a given event and 
the resulting distributions are summed over all events. 

Specifically, A^ is calculated by taking in turn each 

positive pion in an event and incrementing a histogram 
of Ay with respect to all the negative pions in that event. 
A^ is then summed over all events. A similar proce- 
dure is followed for A++, A__, and A ^. Eq. [1] is 

then used to calculate B where A^+(-) is the number of 
positive (negative) pions integrated over all events. The 
balance function is calculated for all events in a given 
centrality bin. In the case of non-identified charged par- 
ticle pairs, relative pseudorapidity {Ar]) is used. Balance 
functions using other variables are presented including 
the relative azimuthal angle, A(/), and the Lorentz in- 
variant momentum difference between the two particles, 
ginv Balance functions in terms of Acf) are sensitive to 
flow and jet effects (See Section IIIIPp . Balance func- 
tions in terms of Qinv are sensitive more directly to the 
temperature of the emitting system (See Section UlIBp . 
In addition, balance functions are presented in terms of 
the components of Qinv in the rest frame of the particle 
pair: g'long, in the beam direction; Qout^ in the direction of 
the transverse momentum of the particle pair; and ^side, 
the direction perpendicular to giong and ^out- Note that 
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The width of the balance function is quantified in sev- 
eral ways. For balance functions in terms of A77, A?/, 
and A0, the widths are calculated in terms of a weighted 
average. For example the width of B{Ar]) is calculated 
as 



(A7?) = 
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where B (Arji) is the value of the balance function for 
the relative pseudorapidity bin Ar]i and the sums are 
carried out from a beginning relative pseudorapidity bin 
Slower to an ending bin iupper- The lower bin is chosen to 
minimize contributions from background and final state 
interactions and the upper bin is the highest bin in A77. 
For balance functions in terms of (^inv, we extract the 
width by fitting to a thermal distribution over a range 



in ginv Widths extracted from the measured balance 
functions are presented in Section V. 

The data used in this analysis were measured using 
the Solenoidal Tracker at RHIC (STAR) J37|, l38|. The 
Au+Au data were acquired during Run 7 at RHIC. The 
p+p data were taken during Run 2 and the d+Au data 
were taken during Run 3. The main detector was the 
Time Projection Chamber (TPC) located in a solenoidal 
magnetic field. The magnetic field magnitude was 0.50 
T. Events were selected according to the distance of their 
event vertex from the center of STAR. Events were ac- 
cepted within 1 cm of the center of STAR in the plane 
perpendicular to the beam direction. Events were ac- 
cepted with vertices within 10 cm of the center of STAR 
in the beam direction for Au+Au and within 15 cm for 
p+p and d+Au collisions. 

Minimum-bias data were used in all cases. Minimum- 
bias triggers for the Au-hAu collisions were defined by the 
coincidence of two Zero Degree Calorimeters (ZDCs) [39| 
located ± 18 m from the center of the interaction region, 
along with an online cut on the Vertex Position Detec- 
tors (VPDs) that restricted accepted events to within 5 
cm of the center of STAR in the beam direction. For 
the Au+Au data set, 28 million events were analyzed. 
For p+p and d+Au collisions, the trigger consisted of 
the two ZDCs combined with the Central Trigger Barrel 
(CTB) [40]. One million events were analyzed for the 
p+p data set and for the d+Au data set. For Au-hAu 
collisions, centrality bins were determined using the mea- 
sured charged hadron multiplicity within the pseudora- 
pidity range |7^| < 0.5 as measured in the TPC. The cen- 
trality bins were calculated as a fraction of this multi- 
plicity distribution starting with the highest multiplici- 
ties. The ranges used were 0-5% (most central), 5-10%, 
10-20%, 20-30%, 30-40%, 40-50%, 50-60%, 60-70%, and 
70-80% (most peripheral). For d+Au, three centrality 
bins were used, 0-20%, 20-60%, and 60-100% determined 
by the multiplicity of charged particles originating from 
the primary collision vertex in the Forward Time Projec- 
tion Chamber (FTPC), in the direction of the deuteron 
beam [4l|. Note that the pseudorapidity distribution for 
d+Au is not symmetric around r] = 0. Each centrality 
was associated with a number of participating nucleons, 
^part, using a Glauber Monte Carlo calculation J4ll, |42|. 
For p+p collisions, all multiplicities were used. 

All tracks were required to have a distance of closest 
approach (DCA) to the measured event vertex of less 
than 3 cm. Only charged particle tracks having more 
than 15 measured space points along the trajectory were 
accepted. The ratio of the numbers of reconstructed 
space points to possible space points along the track was 
required to be greater than 0.52. Charged pions and 
charged kaons were identified using the specific energy 
loss, dE/dx, along the track and the momentum, p, of 
the track. Particle identification was accomplished by se- 
lecting particles whose specific energy losses were within 
two standard deviations of the energy- loss predictions for 
a given particle type and momentum. Particle identifica- 



tion for pions (kaons) also included a condition that the 
specific energy loss should be more than two standard de- 
viations away from the loss predicted for a kaon (pion). 
In addition, electrons were excluded from the analysis for 
all cases by requiring that the specific energy loss for each 
track was more than one standard deviation away from 
the energy-loss predictions for electrons. 

We estimated the systematic errors by comparing the 
results from Run 4 at RHIC with the results from Run 
7 at RHIC, in which new tracking software was imple- 
mented. We assign a 5% systematic error on the ex- 
tracted widths for the balance functions in terms of At^ 
and /S.y and a 10% systematic error on the extracted 
widths for the balance functions in terms of g'lnv and A(/). 

In this paper, we present an overview of the acceptance 
and efficiency of STAR in Section |TT1 because the balance 
functions we present here are not corrected for accep- 
tance and efficiency. This section includes detailed track 
cut speciff cat ions. We then present the balance functions 
for all the measured systems in Section Hill We compare 
some of the results with blast-wave model [30] and HI- 
JING (version 1.38) [43] predictions in Section |IVl We 
then extract the widths of the balance functions and ex- 
amine the systematics of these widths in Section |Vl Our 
conclusions are presented in Section [Vll 



II. DATA ACCEPTANCE AND EFFICIENCY 

Here we outline the major acceptance and efficiency 
corrections necessary to compare any model calculation 
with the balance function results presented in this paper. 
The pseudorapidity cut for all cases is \ri\ < 1.0. The po- 
sition of the vertex for each event along the beam direc- 
tion affects the pseudorapidity acceptance of STAR. The 
distribution of event vertices along the beam direction is 
shown in Fig.[Tl The solid curve in Fig. [1] corresponds to 
a Gaussian fit with a mean of -0.27 cm and a standard 
deviation of 6.81 cm. The distributions of event vertices 
in the beam direction for p+p and d+Au have a standard 
deviation of approximately 25 cm. 

For the balance functions for all charged particles, we 
used a pt cut of 0.2 < pt < 2.0 GeV/c. For identified 
particles, we used a pt cut of 0.2 < pt < 0-6 GeV/c. For 
the high pt measurements for B{A(j))^ we used a pt cut 
of 1.0 < Pt < 10.0 GeV/c. The DC A cut of 3 cm par- 
tially suppressed particles resulting from weak decays. 
The probability of accepting a charged particle in the 
fiducial volume of the TPC (including particle decay) is 
90% for charged particles with pt > 0.2 GeV/c. The 
efficiency for reconstructing a charged pion in our ac- 
ceptance varies from 80% in central collisions to 95% in 
peripheral collisions. More details can be found in Refs. 
Ell, Q? ^^^ O- We also suppressed electrons, result- 
ing in the removal (< 5%) of pions in the momentum 
range 0.20 < p < 0.25 GeV/c. The electron cuts re- 
moved approximately 30% of the identiffed kaons in the 
momentum range 0.4 < p < 0.8 GeV/c. To check these 
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FIG. 1. (Color online) The distribution of the reconstructed 
position of the event vertex along the beam direction for 
events from Au+Au collisions at ^snn = 200 GeV. The solid 
curve is a Gaussian fit with a mean of -0.27 cm and a standard 
deviation of 6.81 cm. 



acceptance and efficiency corrections, we present balance 
functions based on 90k central HIJING events passing 
our event cuts that have been passed through GEANT 
and full event reconstruction. We compare those results 
with our ffltered HIJING calculations in Fig. [51 Filtered 
means that we apply our acceptance cuts in r] and pt 
as well as the efficiency cuts listed above. In addition, 
we present the filtered HIJING calculations with no ef- 
ficiency correction (e = 1), but with all acceptance cuts 
applied. We see that the filtered HIJING results are sim- 
ilar to the full GEANT-filtered HIJING results within 
errors. The widths of all three sets of HIJING data are 
the same within errors. 



III. BALANCE FUNCTIONS 

The balance functions B{Ar]) and B{Ay) can be 
related to the correlation in rapidity of produced 
charge/anti-charge pairs. By comparing PYTHIA cal- 
culations for p+p collisions with the results of a model 
describing a pion gas in which the opposite-charge pion 
pairs are assumed to be created together in space-time, 
the authors of Ref. [30] show that the balance functions 
from p+p collisions were wider than those from a ther- 
mal model. Furthermore, they show that the transport 
model RQMD (Relativistic Quantum Molecular Dynam- 
ics) [46| , in which the hadrons are created during the first 
1 fm/c after the collision, predicts that the balance func- 
tion is wider in more central Au+Au collisions, which is 
the opposite of the experimental trend [23| . We further 
observe that the transport model UrQMD [3i| predicts 
that the widths of the balance function in terms of Arj 
and Ay shows no centrality dependence for s/snn = 200 
GeV Au+Au collisions (see Section V). 

The authors of Ref. [30| make the point that the ob- 
served width of the balance function in terms of relative 
rapidity, a^ , is a combination of the rapidity spread in- 
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FIG. 2. (Color online) Calculated balance functions for all 
charged particles from central Au+Au collisions at ^/sNN 
= 200 GeV using HIJING. The open circles depict HIJING 
events passed through GEANT and event reconstruction. The 
open squares show HIJING events filtered with the acceptance 
and efficiency cuts described in the text. The open triangles 
show HIJING events filtered with the acceptance cuts only. 
When not shown, the statistical errors are smaller than the 
symbol size. 



duced by thermal effects, crtherm, and the separation of 
the balancing partners of the charge/anti-charge pair in 
coordinate space. The authors of Ref. [27] stated this 
relationship as <j^ = cr^herm + 4/31n(r/ro), where /3 is a 
diffusion constant, r is the proper time after the initial 
collision when the charge/anti-charge pair is created, and 
To is a characteristic time on the order of 1 fm/c. After 
the initial collision, the width of the balance function de- 
creases because the thermal width narrows as a result of 
cooling, while diffusion tends to increase the width of the 
balance function. If production of the charge/anti-charge 
pairs occurs at early times, then scattering and expansion 
affects the partners of the charge/anti-charge pair during 
the entire lifetime of the system. The diffusion term is 
then large and significantly broadens the observed bal- 
ance function. If the production of charge/anti-charge 
pairs occurs late, the time during which the partners of 
the charge/anti-charge pair are exposed to scattering and 
expansion is small, which makes the effect of diffusion 
negligible. Thus, in the case of late production of the 
charge/anti-charge pairs, the width of the balance func- 
tion is determined by the thermal width. In Ref. [47], 
the dependence of these model calculations on delayed 
hadronization is demonstrated for a range of model as- 
sumptions. The model calculations show that the longer 
hadronization is delayed, the narrower is the balance 
function. 

In this section, we show the measured balance func- 
tions for p+p, d+Au, and Au+Au collisions at -s/snn = 
200 GeV. We present balance functions for all charged 



particles, for charged pions, and for charged kaons. 
Throughout this paper, plotted balance functions based 
on Eq. [T]for data have been corrected by subtracting the 
balance functions calculated using mixed events. This 
subtraction corrects for differences between the accep- 
tances for positive and negative particles. Mixed events 
are created for each colliding system by grouping the 
events according to bins in centrality and bins in the po- 
sition of the reconstructed vertex of the event along the 
beam direction. For the Au+Au data set, ten centrality 
bins and five vertex bins were used. For the p+p data, 
five bins in event vertex position were used. No mixed 
events were created for the d+Au results because we only 
present results for B{Ar]) for all charged particles, which 
did not require mixed event subtraction. 

A set of mixed events is created by taking one track 
from an event, selected according to the bin in centrality 
and the bin in event vertex position. A mixed event in- 
cludes no more than one track from any observed event. 
This new mixed-event data set has the same number of 
events with the same multiplicity distribution as the orig- 
inal data set but all correlations are removed. B{Ar]) and 
B{Ay) calculated from mixed events are always zero for 
all centralities and for all A77 and Ay. However, the bal- 
ance functions in terms of A^ calculated using mixed 
events are not always zero. The difference between the 
the behavior of positively charged particles and nega- 
tively charged particles crossing the boundary between 
TPC sectors causes the balance functions in terms of A^ 
calculated with mixed events to be non-zero. This effect 
is most pronounced in central collisions where the particle 
density is the highest. These variations of B{A(j)) corre- 
spond to multiples of the azimuthal separation of the sec- 
tor boundaries of the TPC (A^ = 27r/12 = 0.52). Resid- 
ual effects can still be seen in balance functions in terms 
of A^ in the most central bins even after mixed event 
subtraction at A(j) values corresponding to multiples of 
the azimuthal separation of the TPC sector boundaries. 

For most of the measured systems, we also present bal- 
ance functions calculated from shuffled events. These 
shuffled events are produced by randomly shuffling the 
charges of the particles in each event. The shuffled events 
thus have all the momentum correlations and the total 
charge observed in the original event, but the charge- 
momentum correlations are removed. Because shuf- 
fling uniformly distributes a particle's balancing partner 
across the measured phase space, B{Ar]) and B{Ay) cal- 
culated using shuffled events can be used to gauge the 
widest balance functions that one can measure using the 
STAR acceptance for the system under consideration. 
Balance functions calculated with shuffled events have 
the same integral as the balance functions calculated with 
the original events. One exception for the shuffled events 
relates to balance functions calculated using low mul- 
tiplicity events, speciflcally the results for B{Ay) and 
B{qinv) for charged kaon pairs. The balance functions 
calculated by shuffling low multiplicity events are not 
signiflcantly different from the original events, because 
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FIG. 3. (Color online) The balance function in terms of Ar/ 
for all charged particle pairs from Au+Au collisions at ^/sNN 
— 200 GeV for nine centrality bins. 
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FIG. 4. (Color online) The balance function in terms of A^ 
for identified charged pion pairs from Au+Au collisions at 
yJsNN — 200 GeV for nine centrality bins. 



exchanging the positive and negative balancing partners 
has no effect on the resulting balance function. There- 
fore, in the case of low multiplicity events, we create the 
shuffled events by sampling the parent distributions for 
the variable in question. The resulting shuffled balance 
function using sampling has an integral equal to one. The 
shuffled balance functions using sampling are scaled by 
the integral of the original balance function. We veri- 
fied that the shuffled events created using the sampling 
technique agree with the shuffled data in the case of high 
multiplicity events, specifically for B{Ay) and ^(g'inv) for 
charged pion pairs. 



A. Balance Functions in Terms of Ar] and Ay 

1. Au+Au at ^/sj^ = 200 GeV 

Fig. [3] shows the balance function in terms of Arj for 
all charged particles from Au+Au collisions at -s/snn = 
200 GeV for nine centrality bins. The balance function 
gets narrower as the collisions become more central. The 
balance function for mixed events is zero for all centrali- 
ties and A77. The balance function for shuffled events is 
significantly wider than the measured balance functions. 
Model predictions show that inter-pair correlations (e.g. 
HBT and final state interactions) should be significant 
for At/ <0.1 [29]. 

Figs, m and [5] show the balance functions for identi- 
fied charged pion pairs and kaons pairs, respectively, for 
Au+Au collisions at y^iiviv = 200 GeV for nine centrality 
bins as a function of the relative rapidity. The balance 
function for identified pion pairs gets narrower in central 
collisions. The lower magnitude of the balance function 
for pion pairs and kaon pairs compared with the balance 
function for all charged particles is due to the fact that 
the efficiency of observing an identified pion or a kaon is 
lower than for unidentified charged particles. The bal- 
ance function calculated from mixed events is zero for all 
centralities and Ay for both pions and kaons. The bal- 
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FIG. 5. (Color online) The balance function in terms of Ay 
for identified charged kaon pairs from Au+Au collisions at 
^snn = 200 GeV for nine centrality bins. 



ance functions calculated using shuffled events are sub- 
stantially wider than the measured balance functions. 
The discontinuity in B{Ay) for kaons around Ay = 0.4 
visible at all centralities is due to decay, which was 
verified using HIJING calculations. Model predictions 
show that inter-pair correlations should be significant for 
Ay < 0.2 [29|. These effects scale with the multiplicity 
and thus are more apparent in central collisions. 



2. p+p and d+Au at ^/^sWn = 200 GeV 

To investigate the system-size dependence of the bal- 
ance function and to provide a nucleon-nucleon reference 
for the balance functions extracted from Au+Au colli- 
sions, we measured the balance functions for p+p and 
d+Au collisions at ^/sNN = 200 GeV. Fig. [6] shows the 
balance functions for all charged particles for p+p colli- 
sions at ^/s = 200 GeV. The balance functions for p+p 
collisions are integrated over all observed event multiplic- 
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FIG. 6. (Color online) The balance function for p+p collisions 
at v^ = 200 GeV. The top panel shows the balance function 
for all charged particles in terms of Ar]. The bottom panel 
gives the balance function for charged pion pairs and charged 
kaon pairs in terms of A^. 



ities to allow comparison with centrality-selected d+Au 
and Au+Au results. Note that the width of the balance 
function in terms of At^ for p+p collisions is independent 
of the multiplicity of tracks in the event. The top panel 
of Fig. [6] shows the balance function for all charged par- 
ticles in terms of Arj. In the bottom panel of Fig. [6l the 
balance functions are shown for identified charged pion 
pairs and identified charged kaon pairs in terms of Ay 
from p+p collisions at ^/s = 200 GeV. The balance func- 
tion for mixed events is zero for all Arj and all Ay. The 
observed shapes of the balance functions for the identified 
charged pions and kaons are similar to those observed in 
peripheral (70 - 80%) Au+Au collisions. The fact that 
the balance function for kaon pairs has a lower magnitude 
than the balance function for pion pairs reflects the lower 
efficiency for identifying charged kaons versus identifying 
charged pions in STAR. 

Fig. [3 shows the balance functions in terms of A77 for 
all charged particles from d+Au collisions at -s/snn = 
200 GeV for three centrality bins, 0-20%, 20-60%, and 
60-100%. 



B. Balance Functions in Terms of ^inv 

The balance function in terms of Ar] and Ay is ob- 
served to narrow in central collisions and model calcula- 
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FIG. 7. (Color online) The balance function in terms of Ar/ 
for all charged particles from d+Au collisions at ^/sWN — 200 
GeV for three centrality bins. 



tions have been used to interpret this narrowing in terms 
of delayed hadronization [27-30]. However, in a thermal 
model, the width of the balance function in terms of Ar] 
and Ay can be influenced by radial flow. In the absence 
of detector efficiency and acceptance considerations, the 
width of the balance function in terms of the Lorentz in- 
variant momentum difference between the two particles, 
g'inv, is determined solely by the breakup temperature, if 
the balancing charges are emitted from the same position 
in coordinate space. However, when detector acceptance 
is taken into account, some dependence on collective flow 
is introduced ^]. Thus, analyzing the balance function 
in terms of ginv avoids some of the complications associ- 
ated with collective flow, and the balance function calcu- 
lated with a breakup temperature should be the narrow- 
est possible balance function if the particles are assumed 
to be emitted from the same position in coordinate space. 
In addition, contributions to the balance function from 
the decay of particles are more identifiable. For exam- 
ple, the decay of K'g produces a sharp peak in B{qinv) 
for charged pions, while the contribution to B{Ay) for 
charged pions from the decay of Kg is spread out over 
several bins in Ay. 

To study balance functions in terms of qinv, we use 
identified charged pions and identified charged kaons. For 
pion pairs, we observe a peak from the decay K^ -^ 
7r+ + 7r~. For kaon pairs, we observe a peak from the 
decay (j) -^ K^ + K~ . These peaks are superimposed 
on the balance function of correlated charge/anti-charge 
pairs not resulting from the decay of a particle. 



L Au+Au at ^/sJ7N = 200 GeV 

Fig.[8]shows the balance function for identified charged 
pions in terms of qi^v for Au+Au collisions at -s/snn = 
200 GeV for nine centrality bins. These balance functions 
have been corrected by subtracting the balance functions 
calculated using mixed events. These mixed events are 
not zero for all qi^v because of differences in the track- 
ing at TPC sector boundaries for opposite charges. The 
balance functions calculated for mixed events integrate 
to zero as one would expect and the subtraction of the 
mixed events from the measured balance functions does 
not affect the integral of the resulting balance functions. 
At each centrality, a peak is observed corresponding to 
charged pion pairs resulting from Kg ^ 7r+ + 7r~ . The 
sohd curves represent a fit consisting of the sum of two 
terms. The first term consists of a non-relativistic ther- 
mal distribution of the form 



B{qinv) = aqf^v^' 



.7(2^2) 



(3) 



where a is a constant, the pre- factor qf^^ accounts for 
the phase-space effect, and a is a width parameter. The 
second term of the fit is a Gaussian distribution in qi^v 
describing the K^ decay. Note that no peak from the 
decay of the p^ is visible in central collisions around ginv 
= 0.718 GeV/c where one would expect to observe the 
p^. This non-observation of the p^ is in contrast to HI- 
JING, which predicts a large p^ peak, as is demonstrated 
in Section |IVl The p^ peak is visible in the most pe- 
ripheral collisions, which is consistent with our previous 
study of p^ production at higher pt [48]. The authors 
of Ref. [30| attribute the apparent disappearance of the 
p^ in central collisions to the cooling of the system as 
it expands, which lowers the production rate of p^ com- 
pared with pions. The measured balance functions for 
pions are distinctly different from the balance functions 
calculated using shuffled events. In particular, the sharp 
peak from the K^ decay is not present in the balance 
functions calculated using shuffled events. 

HBT/Coulomb effects are visible for ginv < 0.2 GeV/c 
in Fig. [HI Fig. [9] shows the balance function over the 
range of < ginv < 0.2 GeV/c for the most central bin 
(0 - 5%) and the most peripheral bin (70 -80%). The 
Coulomb force pulls opposite charges closer together and 
pushes same charges apart, leading to an enhancement 
of opposite-sign and a suppression of same-sign pairs at 
small ginv This effect leads to a rise in the balance func- 
tion at small ginv, which is larger in central collisions, 
where the long-range Coulomb force affects more parti- 
cles [30]. In peripheral collisions, because the Coulomb 
interaction is less important and the HBT correction is 
larger because of the smaller source size, the Coulomb en- 
hancement disappears and the balance function becomes 
negative at small ginv [30]. 

Fig. [To] shows the balance function for identified 
charged kaons in terms of q^^^ for Au-hAu collisions at 
^/sNN = 200 GeV in nine centrality bins. These balance 
functions were corrected by subtracting mixed events as 
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FIG. 8. (Color online) The balance function in terms of ^inv 
for charged pion pairs from Au+Au collisions at ^/snn — 200 
GeV in nine centrality bins. Curves correspond to a thermal 
distribution (Eq. [3]) plus K% decay. 
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FIG. 9. (Color online) The balance function in terms of ^inv 
for charged-pion pairs in two centrality bins over the range 
< ^inv < 0.2 GeV/c. 



was done for the charged pion results. At each centrality, 
a peak is observed corresponding to charged kaon pairs 
resulting from (j) -^ K~^ -\-K~ . The solid curves represent 
fits consisting of a non-relativistic thermal distribution 
(Eq. [3]) plus a Gaussian distribution in q^nv for the (f) 
decay. HBT/Coulomb effects at low ginv for kaon pairs 
are not as strong as those observed for pion pairs. The 
measured balance functions are distinct from the balance 
functions calculated from shuffled events. 

Several differences between 5(ginv) for charged pi- 
ons and charged kaons are evident. The observed 
HBT/Coulomb effects at low qinv are much stronger for 
pions than for kaons. The HBT/Coulomb effects for pi- 
ons change dramatically with centrality while the HBT 
effects for kaons are small and change little with central- 
ity. The overall normalization for kaons is lower than the 
overall normalization for pions, reflecting the lower effi- 
ciency for detecting identified kaons. The contribution to 
B{qinv) for pions from K^ decay is approximately 7%, in- 
dependent of centrality. The contribution to B{qinv) for 
kaons from (j) decay is approximately 50%, independent 
of centrality. 
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FIG. 10. (Color online) The balance function in terms of ^inv 
for charged kaon pairs from Au+Au collisions at y/sNN = 200 
GeV in nine centrality bins. Curves correspond to a thermal 
(Eq. [3)) distribution plus cj) decay. 



2. p+p at ^/s = 200 GeV 

Fig. [TT] shows the balance functions in terms of ginv for 
p+p collisions at ^/s = 200 GeV. Fig. [TTk shows the bal- 
ance function for charged pion pairs and Fig. [TTb shows 
the balance function for charged kaon pairs. The solid 
curves are thermal fits (Eq. [3|) plus a peak for K^g and 
p^ decay in the case of charged pions, and for (j) decay in 
the case of charged kaons. The thermal fit does not re- 
produce the charged pion results, while it works well for 
the charged kaon data. The mass of the p^ used in the fit 
shown for pion pairs was assumed to be 0.77 GeV/c^. A 
better fit can be obtained if the mass of the p^ is lowered 
by 0.04 GeV/c^, as was observed previously in studies of 
p^ production in p+p collisions at ^/s = 200 GeV |48j . 
This fit is shown as a dashed curve in the upper panel of 
Fig.[TTJ Note that the p^ peak visible in B{qinv) for pions 
from p+p collisions is not observed in B{qi^y) for pions 
from central Au+Au collisions, but is observed for pions 
from peripheral Au+Au collisions, as shown in Fig. [51 



C. Balance Function in Terms of Components of 

^inv 

Here we present results for the three components of 
Qinv These components are q\ong^ the component along 
the beam direction; g'out, the component in the direction 
of the transverse momentum of the observed pair; and 
^side, the component perpendicular to giong and ^out- 

Analysis of the balance function for these three com- 
ponents can address the question of what causes the bal- 
ance function to narrow in central Au+Au collisions. In a 
thermal model where the balancing particles are emitted 
from the same position in coordinate space, the widths 
would be identical for the three components. On the 
other hand, charge separation associated with string dy- 
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FIG. 11. (Color online) The balance function in terms of 
qinv for charged pion pairs [part a)] and charged kaon pairs 
[part b)] from p+p collisions at ^/s — 200 GeV integrated 
over all multiplicities. Solid curves correspond to a thermal 
distribution (Eq. [3|) plus Kg and p^ decay for pions and (j) 
decay for kaons. The dashed curve for pions represents a fit 
to a thermal distribution (Eq. [3|) plus K^ decay and p° decay, 
with the y9° mass shifted down by 0.04 GeV/c^. 



namics should result in balance functions that are wider 
in giong than in gside or ^out |2y, [30| . Also because the 
velocity gradient is much higher in the longitudinal direc- 
tion, diffusion should broaden the balance function more 
in (?iong H- 

Figs. [121 131 and [TH show the balance functions for 
charged pion pairs from Au+Au collisions at -^/snn = 
200 GeV in terms of q\ong^ Qout^ and q^ide respectively. 
The balance functions calculated using mixed events are 
subtracted from the measured balance functions. The 
balance functions for all three components are narrower 
in central collisions than in peripheral collisions. 

The balance functions in terms of q^ide do not look 
like those measured using giong or g'out because the lower 
momentum cut-off of STAR strongly affects 5(^side) for 
^side < 0.38 GeV/c, which underscores the importance of 
performing comparisons with models that have been put 
through detailed efficiency and acceptance filters. 



D. Balance Functions in Terms of Acf) 

The balance function in terms of A0 may yield in- 
formation related to transverse flow at freeze-out l49| 
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FIG. 12. (Color online) The balance function in terms of ^long 
for charged pion pairs from Au+Au collisions at ^/sNN — 200 
GeV in nine centrality bins. 
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FIG. 13. (Color online) The balance function in terms of ^out 
for charged pion pairs from Au+Au collisions at ^/sNN — 200 
GeV in nine centrality bins. 
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FIG. 14. (Color online) The balance function in terms of ^side 
for charged pion pairs from Au+Au collisions at ^ysNN = 200 
GeV for nine centrality bins. 
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FIG. 15. (Color online) The balance function in terms of 
A(/) for all charged particles with 0.2 < pt < 2.0 GeV/c from 
Au+Au collisions at ^/snn = 200 GeV in nine centrality bins. 
The closed circles represent the real data minus the mixed 
events. 



and may be sensitive to jet production. One might ex- 
pect that jet-like phenomena would involve the emis- 
sion of correlated charge/anti-charge pairs at small rela- 
tive azimuthal angles. We present balance functions for 
all charged particles with 0.2 < pt < 2.0 GeV/c from 
Au+Au collisions at ^/sNN = 200 GeV as a function of 
the relative azimuthal angle, A^. In addition, we present 
B{A(j)) for all charged particles with 1.0 < pt < 10.0 
GeV /c to enhance any possible jet-like contributions to 
the balance function. 

Fig. [15] shows the balance functions as a function of 
A^ for all charged particles with 0.2 < p^ < 2.0 GeV/c 
in nine centrality bins. The balance functions for mixed 
events were subtracted. Note that some structure in A^ 
related to the sector boundaries of the STAR TPC is still 
visible after the subtraction of the mixed events. We ob- 
serve a peaking at A^ = in central collisions, while 
in peripheral collisions, the balance functions are almost 
flat. Fig. [15] also shows the balance functions calculated 
using shuffled events. The balance functions from shuf- 
fled events are constant with A0 and show no centrality 
dependence. 

To augment this result. Fig. [16] presents balance func- 
tions in which we use only particles with 1.0 < pt < 10.0 
GeV/c. For this case, we see that the measured bal- 
ance functions vary little with centrality. Again the bal- 
ance functions calculated with shuffled events are con- 
stant with A0 and show no centrality dependence. HI- 
JING calculations for 5(A^) for all charged particles 
with 0.2 < Pt < 2.0 GeV/c exhibit little dependence 
on A(/), while HIJING calculations for particles with 
1.0 < Pt < 10.0 GeV /c are peaked at A(/) = 0, suggesting 
that the balance functions for this higher p^ range show 
jet-like characteristics. 

The dramatically tight correlations in A^ in central 
collisions of Au+Au shown in Fig. [15] are qualitatively 
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FIG. 16. (Color online) The balance function in terms of 
Acf) for all charged particles with 1.0 < pt < 10-0 GeV/c from 
Au+Au collisions at y/sNN — 200 GeV in nine centrality bins. 
The closed circles represent the real data minus the mixed 
events. 



consistent with the radial flow of a perfect liquid. In a 
liquid with very short mean free path, the balancing par- 
ticles would remain in close proximity throughout the re- 
action. A large mean free path, which would necessitate 
a large viscosity, would damp the correlations in A^ [50] . 
This trend is also consistent with a picture where charges 
are not created until after the flow has been established. 



IV. COMPARISON WITH MODELS 

Fig. [T71 compares the measured balance function 
B{Ay) for charged pion pairs from central collisions of 
Au+Au at y/sNN = 200 GeV to the predictions of the 
blast- wave model [30| and to filtered HIJING calculations 
taking into account acceptance and efficiency. The blast- 
wave model includes radial flow, emission of charge/anti- 
charge pairs of particles close together in space and time, 
resonances, HBT and Coulomb effects, strong force ef- 
fects, inter-domain interactions, and a STAR experimen- 
tal filter. The blast-wave calculations shown in Fig. [T71 
include the acceptance cuts in the current paper. The re- 
sulting absolute predictions of the blast- wave model agree 
well with the measured balance function. In contrast, 
the balance function predicted by HIJING is significantly 
wider than the measured balance function. The widths 
of the balance functions predicted by the blast- wave and 
HIJING are compared with the experimental values in 
Fig. [201 

The width of the balance function predicted by the 
blast-wave model is close to the width observed in cen- 
tral collisions. The blast- wave model assumes that the 
charge/anti-charge pairs of particles are created close to- 
gether in space and at the same time, and contains no 
scattering or longitudinal expansion that would widen 
the balance function in terms of Ay. Thus, the agreement 
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FIG. 17. (Color online) The balance function in terms of Ay 
for charged pions from central collisions of Au+Au at ^/sNN 
= 200 GeV compared with predictions from the blast-wave 
model from Ref. 30] and filtered HIJING calculations taking 
into account acceptance and efficiency. 



of the predicted width from the blast- wave model and the 
data is consistent with the idea of delayed hadronization 
in that delayed hadronization in central collisions would 
minimize the contribution of diffusion effects to the width 
of the balance function. 

The balance function in terms of ginv provides the most 
direct way to study the dependence of the balance func- 
tion on temperature. Fig. [18] compares the balance func- 
tion in terms of Qinv for charged pion pairs from central 
collisions of Au+Au at ^/sNN = 200 GeV to the pre- 
dictions of the blast- wave model and to filtered HIJING 
calculations. For the blast- wave model calculations, HBT 
is not included and the decays of the K'^ and p^ are not 
shown. The solid curve for the data represents a fit com- 
prised of a thermal distribution (Eq. [3]) plus K^ decay. 
The dashed curve for the blast-wave model calculations 
represents a thermal fit (Eq. [3j). The dotted curve for 
the HIJING calculations represents a thermal distribu- 
tion (Eq. [3]) plus p^ decay. All the fits are carried out 
over a range in ginv that is not affected by HBT/ Coulomb 
effects. The width extracted from the thermal fit to 
the blast-wave model calculations is compared with the 
width extracted from experimental data in Fig. [2TJ The 
blast-wave model reproduces the observed width in cen- 
tral collisions. The HIJING calculations show a strong 
p^ peak that is not present in the data. 

Future analyses should be able to disentangle the ef- 
fects of cooling and diffusion in driving the narrowing 
of the balance function. Diffusive effects should largely 
manifest themselves in the q\ong variable because the ini- 
tial velocity is in the longitudinal direction and some cre- 
ation mechanisms, such as strings, preferentially separate 
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FIG. 18. (Color online) The balance function in terms of qh 



for charged pions from central collisions of Au+Au at ^/snn 
= 200 GeV compared with predictions from the blast-wave 
model from Ref. 30] and predictions from filtered HIJING cal- 
culations including acceptance and efficiency. For the blast- 
wave calculations, HBT is not included and the decays of the 
Kg and p^ are not shown. 



the pairs in the longitudinal direction. 
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FIG. 19. (Color online) The balance function width (A?7) for 
all charged particles from Au+Au collisions at ^/sWN = 200 
GeV compared with the widths of balance functions calcu- 
lated using shuffled events. Also shown are the balance func- 
tion widths for p+p and d+Au collisions at ^/sNN = 200 GeV. 
Filtered HIJING calculations are also shown for the widths of 
the balance function from p+p and Au+Au collisions. Fil- 
tered UrQMD calculations are shown for the widths of the 
balance function from Au+Au collisions. 



V. BALANCE FUNCTION WIDTHS 

The balance functions presented in the previous sec- 
tion provide insight into the correlation of charge/anti- 
charge pairs in collisions at RHIC. This approach com- 
plements the approach of studying these phenomena us- 
ing charge-dependent correlation functions in two dimen- 
sions, (A77, A^) [18, 19]. The balance function can be 
related to these correlation functions and to other two- 
particle observables. B{Ay) can be interpreted as the dis- 
tribution of relative rapidities of correlated charge/anti- 
charge pairs. The width of B{Ay) then can be used to 
determine whether correlated charge/anti-charge pairs of 
particles are emitted close together or far apart in rapid- 
ity. The width of the balance function B{qinv) can be 
used to study thermal distributions because this balance 
function can be related to the temperature, and is largely 
unaffected by any radial expansion. 

To quantify the evolution of the balance functions 
B{Ay) and B{Ar]) with centrality, we extract the width, 
{Ay) and (A77), using a weighted average (Eq. [2]). For 
B{Ar])^ the weighted average is calculated for 0.1 < Ar] < 
2.0 and for B{Ay)^ the weighted average is calculated for 
0.2 <Ay< 2.0. 

Fig. [19] shows the balance function widths for all 
charged particles from Au+Au, d+Au, and p+p colli- 
sions at ^/sNN = 200 GeV plotted in terms of the number 
of participating nucleons, A^part- In addition, we present 



the widths of the balance functions from Au+Au col- 
lisions for shuffled events. The widths of the shuffled 
events are considerably larger than those from the mea- 
sured data and represent the largest width we can mea- 
sure using the STAR acceptance for the system under 
consideration. 

The balance function widths scale smoothly from p+p 
through the three centrality bins for d+Au and down to 
the nine Au+Au collision centrality data points. This 
figure also shows filtered HIJING calculations for p+p 
and Au+Au calculations for HIJING and UrQMD. The 
HIJING calculations for p+p reproduce the measured 
width. The Au+Au HIJING and UrQMD calculations, 
however, show little centrality dependence and are com- 
parable to those calculated from the HIJING p+p sim- 
ulations. This is despite the fact that HIJING does not 
predict any appreciable radial flow while UrQMD pre- 
dicts radial flow in Au+Au collisions but less than that 
observed experimentally. This radial flow should produce 
a narrower balance function in central collisions where ra- 
dial flow is the largest, while hadronic scattering should 
lead to a wider balance function. The fact that the mea- 
sured widths from Au+Au collisions narrow in central 
collisions is consistent with trends predicted by models 
incorporating late hadronization [27, 30]. 

Fig. [20] presents the widths of the balance func- 
tion, B{Ay)^ for identified charged pions and identified 
charged kaons from p+p collisions at ^/s = 200 GeV and 
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FIG. 20. (Color online) The balance function widths for iden- 
tified charged pions and charged kaons from Au+Au collisions 
at ^/s]^ = 200 GeV and p+p colhsions at ^ = 200 GeV. Fil- 
tered HIJING calculations are shown for the same systems. 
Filtered UrQMD calculations are shown for Au+Au. Also 
shown is the width of the balance function for pions predicted 
by the blast- wave model of Ref. ]3u\ . 



Au+Au collisions at s/snn = 200 GeV. Also shown are 
filtered HIJING and UrQMD calculations. For charged 
pions, the measured balance function widths for Au+Au 
collisions get smaller in central collisions, while the fil- 
tered HIJING and UrQMD calculations for Au+Au again 
show no centrality dependence. The HIJING calculations 
for p+p collisions reproduce the observed widths. 

In contrast, the widths of the measured balance func- 
tion for charged kaons from Au+Au collisions show little 
centrality dependence. The extracted widths for charged 
kaons are consistent with the predictions from filtered 
HIJING calculations and are consistent with the p+p re- 
sults. The widths for charged kaons predicted by UrQMD 
are somewhat larger than the data. The agreement with 
HIJING and the lack of centrality dependence may in- 
dicate that kaons are produced mainly at the beginning 
of the collision rather than during a later hadronization 
stage [27]. The larger widths predicted by UrQMD for 
kaons may reflect the hadronic scattering incorporated 
in UrQMD, although the statistical errors are large for 
both the data and the model predictions. 

Fig. [21] shows the widths extracted from 5(ginv) for 
identified charged pions and kaons from Au+Au collisions 
at ^/sNN = 200 GeV and p+p collisions at ^/s = 200 
GeV using a thermal distribution (Eq. [3]) where a is the 
width. The widths for the pions are somewhat smaller 
than the widths for the kaons, although the kaon widths 
have a large statistical error. This width is related to the 
temperature of the system when the pions and kaons are 



formed. Filtered HIJING calculations show no centrality 
dependence and predict a difference between the widths 
for pions and kaons. The widths predicted by UrQMD for 
pions are smaller than those predicted by HIJING but are 
still larger than the measured widths. In addition, the 
widths predicted by UrQMD for pions seem to show a 
centrality dependence, although it is not as strong as that 
for the data. The widths predicted by UrQMD for kaons 
show no centrality dependence and agree with HIJING. 

For a thermal system in the non-relativistic limit {m ^ 
T), the balance function has the functional form given in 
Eq. [3] where a = \^2mT. For kinetic freeze-out tem- 
peratures T ~ 0.1 GeV 0, kaons are non-relativistic, 
and this functional form was seen to describe the balance 
function in Fig. [TOl Indeed, as seen in the right panel of 
Fig. EU the evolution in the width of the balance func- 
tion may be understood in terms of the evolution of the 
freeze-out temperature as a function of centrality [45[ . 

In the ultra-relativistic case {m <C T), the balance 
function from a thermal system is exponential rather 
than Gaussian, B{qinv) ^ Q\nw^~^'''^^^ - The proper func- 
tional form for pions, being neither non-relativistic nor 
ultra-relativistic, is more complicated. Indeed, we found 
that neither the Gaussian form nor the exponential form 
fully describe the pion balance function in Fig. [51 Thus, 
to get a feeling for whether the evolution in freeze-out 
temperature can explain the narrowing of the balance 
function for pions, we turn to numerical calculations. 
Calculations in Ref. [29| show a 27% reduction in the 
Gaussian width of B{q\rw) as the temperature is varied 
from 120 to 90 MeV, the temperatures inferred from fits 
to peripheral and central collisions, respectively [45]. As 
seen in Fig. [211 the measured width for peripheral (cen- 
tral) collisions is 0.33 GeV/c (0.27 GeV/c), a 18% reduc- 
tion. Thus, the centrality evolution in freeze-out tem- 
perature may help explain much of the narrowing of the 
balance function in terms of ginv for pions as well as for 
kaons. However, firm conclusions require more complete 
calculations including all detector effects. 

Fig. [221 shows the widths of the balance functions in 
terms of ^long, ^out, and ^gide for charged pion pairs in 
Au+Au collisions at -s/snn = 200 GeV compared with 
the results of filtered UrQMD calculations. These widths 
were extracted by taking the weighted average over the 
<7iong5 Qout^ and ^gide range from 0.0 to 1.3 GeV/c. The 
width (qside) is larger than (g'long) and (g'out) because 
the lower pt threshold of STAR aflFects it more strongly. 
In the most peripheral collisions, the widths (^long) and 
(^out) are comparable to each other. As the colhsions be- 
come more central, both ((^long) and (qout) decrease. The 
change in (giong) is less than the change of (^out) with in- 
creasing centrality. Thus it seems that the two transverse 
widths, (^out), and (^side), decrease in central collisions 
more strongly than the longitudinal width, (giong)- This 
may imply that string dynamics and diffusion due to lon- 
gitudinal expansion may keep (giong) from decreasing as 
much in more central collisions [30|. The decrease in the 
transverse widths is consistent with the decrease in Tkin 
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FIG. 21. (Color online) The balance function width a extracted from B{qin^/) for identified charged pions and kaons from 
Au+Au collisions at y^iiviv = 200 GeV and p+p collisions at y/s = 200 GeV using a thermal fit (Eq. [S]) where a is the width. 
Filtered HIJING and UrQMD calculations are shown for pions and kaons from Au+Au collisions at ^/sNN — 200 GeV. Values 
are shown for \/2rnT^ from Au+Au collisions, where m is the mass of a pion or a kaon, and Tkin is calculated from identified 
particle spectra |45[. The width predicted by the blast- wave model of Ref. [30] is also shown for pions. 
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For the lower pt particles, the balance function narrows 
dramatically in central collisions (large positive values of 
(cos(A0))). The narrow balance functions observed in 
central collisions may be a signature of the flow of a per- 
fect liquid, as discussed above. For the higher pi parti- 
cles, (cos (A0)) in Au+Au collisions shows less centrality 
dependence. 

Fig. [231 also shows UrQMD calculations for (cos (A(/))). 
The predictions for the 0.2 < pt < 2.0 GeV /c data set are 
much lower than the measured values, which is consistent 
with the observation that UrQMD underpredicts radial 
flow. The predictions for (cos(A^)) for the 1.0 < pt < 
10.0 GeV/c data set show no centrality dependence and 
are also much lower than the measured values. 



FIG. 22. (Color online) The widths for the balance func- 
tions for pions in terms of qiong, ^out, and ^side compared with 
UrQMD calculations. 



as the collisions become more central. In the most pe- 
ripheral collisions, the widths predicted by UrQMD are 
consistent with the data. As the collisions become more 
central, the predicted widths decrease slightly, but not 
as much as observed in the data. This is consistent with 
results using the balance function in terms of qinv Addi- 
tional theoretical input is required to draw more conclu- 
sions from the analysis of the balance function in terms 
of the components of ginv 

Fig. [23] shows the weighted average cosine of the rel- 
ative azimuthal angle, (cos(A^)), extracted from the 
balance functions 5(A0) for all charged particles from 
Au+Au collisions at ^/sNN = 200 GeV with 0.2 < pt < 
2.0 GeV/c and 1.0 < Pt < 10.0 GeV/c. The values for 
(cos(A^)) are extracted over the range < A0 < tt. 



VI. CONCLUSIONS 

We have measured balance functions for p+p, d+Au, 
and Au+Au collisions at ^/sNN = 200 GeV for all 
charged particles, identified charged pions, and identi- 
fied charged kaons. We observe that the balance func- 
tions in terms of A77 for all charged particles and in 
terms of Ay and ginv for charged pions narrow in central 
Au+Au collisions. This centrality dependence is consis- 
tent with trends predicted by models incorporating de- 
layed hadronization. The balance functions B{Ar]) and 
B{Ay) can be affected by radial flow while the balance 
function 5(ginv) is largely unaffected by the implied ref- 
erence frame transformation. We observe that the sys- 
tem size dependence of the width of the balance function 
for charged particles scales with A/'part as was observed 
at ./sj^ = 17.3 GeV [34]. In contrast, HIJING and 
UrQMD model calculations for the width of the balance 
function in terms of Ay or At^ show no dependence on 
system size or centrality. 
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FIG. 23. (Color online) The weighted average cosine of the 
relative azimuthal angle, (cos(A(/))), extracted from B(A(j)) 
for all charged particles with 0.2 < pt < 2.0 from Au+Au 
collisions at y/sNN = 200 GeV and from all charged particles 
with 1.0 < Pt < 10.0 GeV/c, compared with predictions using 
filtered UrQMD calculations. 



For charged kaons we observe that the width of the 
balance function B{Ay) shows little dependence on cen- 
trality for Au+Au collisions at ^/sNN = 200 GeV. 
This lack of dependence on centrality may indicate that 
strangeness is created early in the collision rather than 
in a later hadronization stage. However, the fact that 
the balance function for kaons in terms of q^nv narrows 
in central collisions might be explained by the exclusion 
of the (j) decay in the fits to B{qinv)^ while the (j) decays 
are included in B{Ay). 

For both pions and kaons, the width of the balance 
function in Qinv decreases with increasing centrality. This 
narrowing may be driven largely by the evolution of the 
kinetic freeze-out temperature with centrality. This ex- 
planation is strengthened by the observation that the 
widths of the balance functions for pions in terms of 
the two transverse components of g'lnv, ^out and ^side, de- 



crease in central collisions. However, more quantitative 
conclusions require more complete theoretical studies. 

A comparison with a blast-wave model J30J suggests 
that the balance function B{Ay) for pion pairs in cen- 
tral Au+Au collisions at -^snn = 200 GeV is as nar- 
row as one could expect, as the model assumed that the 
balancing charges were perfectly correlated in coordinate 
space at breakup. This correlation might be explained 
either by having the charges created late in the reaction, 
thus denying them the opportunity to separate in coordi- 
nate space, or having them created early, but maintain- 
ing their close proximity through very limited diffusion. 
Whereas the first explanation is motivated by a picture 
of delayed hadronization, the idea of limited diffusion is 
consistent with the matter having a very small viscosity, 
which also requires a small mean free path. Furthermore, 
both these explanations account for the observation that 
the balance function narrows with centrality, since the 
breakup temperature, which determines the width, falls 
with increasing centrality. The additional information 
provided here concerning the decomposition of the bal- 
ance function into g'out, ^side, and g'long may provide the 
basis for a more stringent test of competing theoretical 
pictures. 
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